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Abstract 

Self-similar solutions to the problem of a special relativistic law of 
motion for thin shells of matter are calculated. These solutions rep- 
resent the special relativistic generalization of momentum conservation 
for the thin layer approximation in classical physics. The analytical and 
numerical results are applied to Supernova Remnant 1987 A. 
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1 Introduction 

The study of the relativistic dynamics of thin shells of matter is a current 
subject of investigation when the framework of general relativity is adopted 
PQ |2j [3j H]. Here, we will explore how the framework of special relativity 
can produce a law of motion which can be compared with expansion data 
of a supernova remnant (SNR). From a classical point of view, the temporal 
evolution of the SNR is modelled by adopting different physical approaches 
that are not invariant under relativistic transformations. As an example, a 
classical analytical solution of the Sedov type (5j[6l[7] is equation (10.27) in [8] 

Rit) = ( ^ )1/5 ' (1) 

where p is the density of the surrounding medium which is supposed to be 
constant, E is the energy of the explosion and t is the age of the SNR. In 
this paper, we show that it is possible to deduce the relativistic law of mo- 
tion starting from the conservation of relativistic momentum. The relativistic 
equation of a SNR is deduced as a non-linear relationship between radius and 
time which can be solved numerically. The self-similar relativistic solutions 
are accurate for sufficiently large time after the formation of the SNR. The 
theory is applied to the SNR connected with Supernova (SN) 1987A. 



2 



L. Zaninetti 



2 Classical and relativistic laws of motion 

The thin layer approximation in classical physics assumes that all the swept-up 
gas accumulates with infinite density in a thin shell just after the shock front. 
The conservation of radial momentum requires that, after the initial radius 

^R 3 pV=^R 3 pV , (2) 

where R and V are the radius and velocity of the advancing shock wave, p is 
the density of the ambient medium and Vq is the initial velocity, see [10] . In 
classical physics, the velocity as a function of radius is: 

V = V (§) 3 , (3) 



the law of motion is: 



i 

Z? = /?n(l + 4-^(*-* )V , (4) 



o / 

where t is time and to is the initial time. In classical physics, the velocity as a 
function of time is: 

_ 3 

V = V (l + 4^-(t-t )^) * . (5) 

Equation (j3J) can also be solved with a similar solution of type R = K(t — t ) a , 
k being a constant, and the classical result is: 



R=n{/f3 R Q 3 c(t-t ^ , (6) 

where (3q = — has been introduced in order to make a comparison with the 
relativistic case. 

Newton's law in special relativity is: 

where F is the force, p the relativistic momentum, m the relativistic mass, itlq 
the mass at rest and V the velocity, see equation (7.16) in [IT]. In the case of 
the relativistic expansion of a shell in which all the swept material resides at 
two different points, denoted by radius R and Ro, we have: 
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Figure 1: Velocity as a function of radius when Rq — 1, c = 1, Vq/c = 0.99. 
(relativistic case, equation (jHJ)) (dashed-line ) and (classical case, equation 
(full-line ) 

where (3o=Vq/c, 0—V/c and c is the velocity of light. The velocity of the 
relativistic expanding shell is: 

= ^2 ( 9 ) 

yjRoW + R 6 - R 6 Po 2 

Figure [1] shows the classical and relativistic behaviors of the velocity as a 
function of radius R. 

The previous formula can be expressed in differential form as: 



iVA) +R 6 ~ R 6 /3 2 dR = R d (3 cdt . (10) 

The integral on the lhs of the previous equation can be evaluated with a first 
transformation u = D l a — -L and x = R , 

Jyfl + fiX 6 dx . (11) 
A second transformation y = x 6 changes the integral into: 



lVTT ^dy . (12) 



6 y 5 / 6 

This integral is of the same type as formula 3.194.1 in [12 



x» 1 ^ = ut* 2 Fi(fJ.,v; 1 + //; ~/3u) 



o (1 + /3 x) u IL 
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Figure 2: Plot of the hypergeometric function 2-^1 ( — 1/2, 1/6; 7/6; 2) as a 
function of x (full line) for —10 < x < +1 . 



where 2 F\ (a, 6; c; z) is a regularized hypergeometric function [131 EH US E2]- 
In our case, v = | and /i = ~. The hypergeometric function is defined by the 
following power series expansion: 



2 F 1 (a,b;c;z) = ±^^z n 

where (w) n is the Pochhammer symbol 

(w) n = w(w + 1) . . . (w + n — 1) 



(14) 



(15) 



(a, b, c) is a triplet of real numbers with c not belonging to the set of negative 
integers and z is a real number < 1. 

We are now ready to integrate formula (fTUl ) and the result is a non-linear 
equation, J^nl, in R- 



NL 



i? 2^-1/2, 1/6; 7/6; 



i? 6 (-i + A> s 



-i? 2 iM-l/2,l/6; 7/6; 



-1 + A) 2 
A 2 



) 



-c(t - t )R 3 p = 



(16) 



From a numerical point of view, the hypergeometric function can be evalu- 
ated with the FORTRAN subroutine HYGFX extracted from [16] and a typical 
plot of 2-^1(1/6, 1/2; 7/6; x) is reported in Figure |2j 
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Figure 3: Radius as a function of time when R = 1, (3 = 0.99, c = 1, (asymp- 
totic relativistic case, equation ffTT|) ) (dot-dash-dot-dash line), (relativistic 
case, equation ffTB"]) ) (dashed-line ) and (classical case, equation @) (full-line) 



Once the hypergeometric function is implemented, we can solve the non- 
linear equation ( JT6l) with the FORTRAN subroutine ZRIDDR from [IT] and 
a typical example is shown in Figure fl3]). Once a numerical expression for the 
relativistic radius is obtained, we can easily obtain a velocity-time relationship 
from equation ([9]) , see Figure HJ 

An approximate solution of equation [9] can be found by imposing R(t) = 
k(t — to) a . The self-similar solutions of the relativistic case under the assump- 
tion R 6 (l - (3D > Rlfil are: 



R(t) = V2 ^ 
\ 



i-A) 2 



{t-k 



(17) 



and 



0{t) = l/Ay/2 



-3/4 c -l 



3 Application to SN 1987A 

The SN 1987A exploded in the Large Magellanic Cloud in 1987. The distance 
of this SN is ~ 50 kpc (163050 ly) and a detailed analysis of the distance, D, 
gives D = 51.4 kpc [TS] and D = 50.18 kpc [19J. After 7987 days, the diameter 
of the SNR was 0.77" and it's velocity « 1412 km/s [20][2TJ. In this Section 
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Figure 4: Velocity as a function of time when R = 1, O = 0.99 and c = 1. The 
dashed-line represents the relativistic case, the dot- dash- dot- dash represents 
the asymptotic relativistic case and the full line the classical case. 

we will adopt year (yr) as a time unit and light year (ly) as a length unit; with 
these natural units c = 1. The radius of the SNR after 21.86 yr is 

R = r 77 x DjsoO.61 ly , (19) 

where r 77 is the radius in arcsec divided by 0.77 and D 50 the distance in pc 
divided by 50000. Next, we attempt to evaluate the initial conditions Ro and 
fio after t = 21.86 yr and R = 0.61 ly are given. The approximate self-similar 
relativistic solution for the radius as represented by equation ( TT71) allows a 
relationship to be determined between these two unknown variables, 

A> = i R \ , (20) 
'R* + 16 RrfcH 2 



and Figure [5] reports such correlated initial conditions. Figure [6] shows the 
behavior of the radius as a function of time after setting initial conditions 
given by equation ( 1201) and extracting the data of SNR 1987A from Figure 2 
in [2]. 

A further test can be done by inserting in formula ffTSl for the self-similar 
relativistic velocity the boundary conditions used to deduce the trajectory, i.e. 
Ro = 0.079, /3q = 0.95 and t = 21.86 yr; the theoretical velocity turns out to 
be 2672 against the observed 1412 *f [21J. The proper time of the world 
line is 

r-r = f yfl - /3(tydt . (21) 
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Figure 6: Radius as a function of time when R = 0.079, (3q = 0.95, c = 1 
(asymptotic relativistic case, equation ( 1T7|) ) (dotted line) with the addition of 
the observed radius with relative error extracted from [21] 
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Log 10 (t -t ) 



Figure 7: The time contraction r of the clock on the expanding surface as 
a function of the decimal logarithm in laboratory time when R = 0.079, 
0o = 0.95 and c = 1. 

This means a relativistic time contraction, r, for the clock which follows the 



expansion 



r -t 
t-tn 



(22) 



and Figure [7] shows such a contraction of the framework which sees the expan- 
sion when the velocity is evaluated as a function of time from equation (Q. 



4 Conclusions 

The introduction of a relativistic framework into the equation of a SNR under 
the hypothesis of the thin layer approximation avoids the paradox of an initial 
velocity greater than the velocity of light. The self-similar solutions for radius 
and velocity, respectively eqns. (fT7|) and ( 118]) . are found under the approxima- 
tion R 6 (l -pi) > Rlfil . The application of these new formulae to SNR1987A 
produces acceptable results. The conservation of classical and relativistic mo- 
mentum adopted here does not take into account the momentum carried away 
by photons. 
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